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Abstract 

An experimental test of the electron energy scale linearities of SNO+ and EJ-301 scintillators was carried out using a 
Compton spectrometer with electrons in the energy range 0.09-3 MeV. The linearity of the apparatus was explicitly 
demonstrated. It was found that the response of both types of scintillators with respect to electrons becomes non-linear 
below ~0.4 MeV. An explanation is given in terms of Cherenkov light absorption and re-emission by the scintillators. 
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1. Introduction 

SNO+ is a multi-purpose neutrino experiment whose 
reach extends to the following areas of neutrino physics: 
neutrinoless double beta decay (with Nd-loaded scintil- 
lator), geo-neutrinos, reactor and low-energy solar neu- 
trinos, as well as supernova neutrinos [1]. It is a ~780- 
tonne liquid scintillator (LS) detector currently under 
construction ~2 km underground at the SNOLAB fa- 
cility in Sudbury, Ontario, Canada. The scintillator will 
be contained in a 12-m diameter spherical acrylic ves- 
sel surrounded by ~9500 photomultiplier tubes (PMTs). 
More details on the layout can be found in 

The LS to be used in SNQh- is linear alkylbenzene 
(LAB) with ~2 g/L of PPG (2,5-diphenyloxazole). In 
the double-beta decay phase, the LAB-PPO will be 
loaded with ~0.1 % Nd (by weight). To achieve the 
goals of the experiment, it is imperative to understand 
the properties of this scintillator down to the SNO-H 
energy threshold, which is expected to be ~0.2 MeV 
for low-energy solar neutrino studies. In particular, for 
electrons, a precise knowledge of the dependence of the 
light yield, X, on their kinetic energy, Ee, is required. 
For scintillators, X(£f ) is commonly assumed to be lin- 
ear, except at very low energies (typically below ~50 
keV), where the effects of ionization quenching appear 
to be important. The decrease in light yield due to ion- 
ization quenching is well-described for all particles by 
Birks' law fl. 
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29 In the last decade, the KamLAND experiment ob- 

30 served, using y-ray calibration sources, that the elec- 

31 tron energy scale of a dodecane-based scintillator was 

32 non-linear and could not be solely described by Birks' 

33 law t4i]. A laboratory investigation of the KamLAND 

34 LS electron energy scale with a Compton spectrome- 

35 ter confirmed these results [5]. It was hypothesized that 

36 Cherenkov light, being absorbed and re-emitted by the 

37 LS, acted as another scintillation light 'source' besides 

38 ionization. This led to a rise in dX/d/s^ in the energy 

39 region above the Cherenkov threshold for electrons. 

40 The primary aim of this work is to test the linearity 

41 of £XEe) for both LAB-PPO and Nd-doped LAB-PPO, 

42 in an energy range relevant to the various goals of the 

43 SNO+ experiment, i.e. up to 3.5 MeV. We verify if the 

44 SNO+ LS displays the same behaviour as KamLAND's, 

45 and quantify the extent of the Cherenkov contribution 

46 to non-linearity. We also test the linearity of EJ-301 

47 scintillator [6], which is equivalent to the commonly- 

48 used NE-213. Details of our apparatus and method are 

49 given in ^ |3] describes a measurement of the non- 
50 linearity of our set-up. The results of our scintillator 

51 linearity tests are presented and discussed in ^ We 

52 summarize and conclude in ^ 



53 2. Experimental method 

54 To measure the scintillation light output as a func- 

55 tion of Ee, we used the Compton scattering of y-rays 

56 of known energies, Ey, to produce mono-energetic elec- 
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Figure 1: Set-up of the appai'atus. 



57 trons within the LS. Ee is given by: 
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where E'^ is the scattered y energy. E'^ is calculated us- 
ing the Compton formula: 
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where is the electron mass, and 6 is the j scatter- 
ing angle (defined in the lower part of fig. [!}. The scat- 
tered J was recorded with a Nal detector placed at angle 
Q. Only events that resulted in coincident pulses from 
the LS and Nal detectors were considered. Two y-ray 
sources of 0.662 and 4.44 MeV allowed values of E^ in 
the range 0.09-3 MeV to be probed with high granular- 
ity. 

Fig . [Bis a photograph of the spectrometer set-up, with 
the LS cell located L5 m from the y source, and a 2" 
Nal detector placed on a trolley with a 0.5-m rotation 
radius. The y-rays were collimated by a 5-mm diameter 
hole drilled through 20 cm of lead shielding. The 0.662 
MeV ys originated from a '^'^Cs source that was small 
enough to fit into the collimator hole. The 4.44 MeV 
source is described below. 
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2.1. A 4.44 MeV y -ray source 

We constructed a high intensity 4.44 MeV y source 
using the '^C(p,p')'^C reaction. A ~2 ;uA 5.7 MeV pro- 
ton beam was supplied by the CENPA tandem Van de 
Graafi" accelerator. Our source chamber, seen in fig. [1] 
was made of aluminium. It consisted of a flange reducer 
section that interfaced to the beam pipe, and an insu- 
lated target holder section connected to a pico-ammeter 
for flux monitoring. A small tantalum aperture located 
at the beam pipe interface ensured that most ys were 
produced at the centre of the carbon target, which was a 
5-mm thick disk of natural carbon. Cooling for the lat- 
ter was not required. We carefully aligned the LS cell, 
y-ray collimator, and the proton beam coUimating and 
focussing apertures with a surveying telescope located 
far away in the 0° direction. 

The '^C(p,p')'^C cross-section was previously mea- 
sured by Barnard et al. and Dyer et al. iJHl- The beam 
was operated at 5.7 MeV to select the first resonance 
peak and reduce backscattered y backgrounds in the ex- 
perimental hall as much as possible. At 2 yuA of 5.7 
MeV proton beam current, and assuming a cross-section 
of 200 mb, we estimated the overall source strength to 
be a few tens of mCi. 

A Nal detector placed along the beamline and 2 m 
from the Carbon target, confirmed the intense 4.44 MeV 
gamma event rate. Two other photons, with energies of 
3.11 and 0.511 MeV, were emitted from '^C(p,p')'^^C, 
and from the decay of '^'N produced in the '^C(p,n)'^N 
reaction. The relative contributions of the 4.44, 3.11 
and 0.511 MeV ys are 89.1, 0.9 and 10%, respectively 
(0. At any angle 0, the values of Eg for recoil elec- 
trons in the LS from the 0.51 1 MeV and 4.44 MeV ys 
were very dififerent, while the 3.11 MeV y interfered 
with the electron signal from the 4.44 MeV line at the 
< 1% level. Therefore, the 0.511 and 3.11 MeV fines 
can be neglected in this work. 

2.2. Scintillator cell 

Two scintillator-PMT assemblies were custom- 
designed and built with all materials in direct contact 
with the active volume being compatible with the LAB- 
based and EJ-301 scintillators. The scintillator was cou- 
pled to a Philips Valvo XP-2262 PMT, housed in mu- 
metal shielding, via a pyrex window. The PMT was 
biased to -1400 V. The 2" x 2" cylindrical scintilla- 
tor container was bored out of a solid piece of stainless 
steel. Its walls were 1 mm thick, with mirror-polished 
inner surfaces. 

Both of our LAB-based LS samples contained 3 g/L 
of PPO. Before filling, the LS was de-oxygenated by 
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Figure 2: Schematic of the electronics chain. 

gently bubbling in Ar gas for ~45 minutes IS- The 
cells were cleaned in an ultra-sound bath, before being 
vacuum-filled with the LS. We used tubes and valves 
that were made entirely of teflon. 

2.3. Electronics and data acquisition 

An NIM-based electronics set-up (shown in fig. |2]l 
recorded the charge of coincident PMT pulses from the 
Nal and liquid scintillator detectors. The pulse height 
spectra from the multi-channel analyzer (MCA) were 
read in by ORCA, a data-acquisition software developed 
at the Universities of Washington and North Carolina 

2.4. Data-taking procedure 

The electronics and PMT were allowed to stabilize 
overnight before data-taking the next day. To avoid 
making corrections due to gain drifts, measurements 
with the 4.44 MeV and 0.662 MeV ys were done on the 
same day, with the '^^Cs runs starting soon after switch- 
ing off the proton beam. For consistency, the LAB- 
PPO and Nd-doped LAB-PPO measurements were per- 
formed using the same scintillator cell and the same 
PMT, which underwent the linearity tests described in 
^ The electronics settings were fixed throughout the 
whole experiment. 

Fig.[3]shows the LS pulse height spectra at seven an- 
gles between 28.5° and 141°. This data was taken with 
the '^^Cs source and LAB-PPO target. At each angle. 
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Figure 3: Top: From left to right, the solid-line histograms are the 
observed electron recoils in LAB-PPO resulting from the 28.5°, 36°, 
46°, 56°, 76°, 96° and 141° scattering of y-rays originating from a 
'^'Cs source. The reconstructed values of Eg at each angle are given. 
The dotted line is the full Compton scattering spectrum. Bottom: The 
LAB-PPO Compton spectra of 0.662, -1.25 and 4.44 MeV ys from 
'"Cs, ^"Co and the '2C(p,p')'-C reaction, respectively. 



155 mono-energetic electron peaks can unambiguously be 

156 located to give X in units of ADC bins. The full Comp- 

157 ton spectrum, drawn in dotted line, is shown for com- 

158 parison. The bottom plot shows two more Compton 

159 spectra, obtained with a ^"Co source and our 4.44 MeV 

160 accelerator source. 



161 3. Linearity of apparatus 

162 The recorded ADC peak positions from the LS can 

163 scale non-linearly with Eg for the following reasons: 

164 (1) the scintillator response is intrinsically non-linear, 

165 (2) the integrated charge from the XP-2262 PMT pulses 

166 did not increase linearly with the number of photons in- 

167 cident on its photocathode, (3) the electronics chain dis- 

168 torted the energy scale, and (4) the apparatus was mis- 

169 aligned, resulting in erroneous values of E^ when using 

170 the Compton equation with scattering angles that were 
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Figure 4: Top: The angulai' dependence of Nal photo-peak positions ™ 

(black markers) of ys from a '^^^Cs source, scattered off a Nd-doped 199 

LAB-PPO target. The dotted Une is a fit to Eq.[3] Bottom: Fit residu- ^oo 
als. 



not corrected for this misalignment. The last three is- 
sues are discussed below. 

3.1. Misalignment and electron energy reconstruction 

The spectrometer misalignment can be studied by fit- 
ting the angular dependence of Nal photo-peak posi- 
tions, /(£■;), from i^^Cs source data to the Compton 
formula, Eq.|2] The fit function was: 



a ■ E-v 



1 + A,(i_cos(0 + 5)) 



+ b 



(3) 



where a and b are constants, and 6 is an angular offset 
to account for the misahgnment. Fig.|4]shows Nal data 
(black markers) fitted to Eq.[3]from 0° to 141°. The un- 
certainty on reading out 6 was estimated to be ±0.25°. 
Except at ^ = 0°, all the data were collected in coinci- 
dence with the LS. When 6 was fixed to 0°, the;\fVNDF 
for the fit was 19.9/15. When it was allowed to float, the 
best fit value for 6 was -0.36° ±0.13°, with an improved 
;\f^/NDF of 7.8/14. To account for this angular off'set, 
we increase our uncertainty in the scattering angle from 
±0.25° to ±0.36°. The residuals for the unconstrained 
fit are shown in the lower plot of fig. |4] 

We performed further checks on the reconstructed 
values using another method, which consists of calibrat- 
ing the Nal energy scale with ^^^Cs, ^*'^Bi and ^^^Eu 
y-ray sources. A fit to a straight line below 0.7 MeV 
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Figure 5: Schematic of the LED driver circuit. 

yielded values of a and b that were consistent with the 
Compton fit. The calibration constants were then used 
to get Ey, and hence E^. This method is independent 
of the one described in |2l because no assumptions are 
made about Compton scattering. The mean percentage 
diff'erence between the values of E^ reconstructed using 
the two methods was -0.4 % for '"Cs data with a LAB- 
PPO target, and 1 % for a Nd-doped LAB-PPO target. 

Note that the second method is applicable only when 
the Nal photo-peaks for the scattered ys are clearly visi- 
ble. For the 4.44 MeV y data, this was not the case, and 
it was assumed Eq.|2] still accurately applied. 

3.2. Measurement of the linearity ofPMTand electron- 
ics chain 

For the LAB-based LS, the linearity of the PMT and 
electronics chain was verified using a system of two se- 
quentially flashing Light Emitting Diodes (LEDs), fol- 
lowing the method described in (l^l- The two LEDs 
(diff'use yellow HLMP-3401) were placed adjacent to 
each other, 10 cm in front of the XP-2662 in a light-tight 
box. The PMT was connected to exactly the same elec- 
tronics (with the same voltage and gain settings) as in 
the scintillator linearity measurements. An LED driver 
(fig. |5]l continuously flashed the LEDs in the following 
sequence consisting of three steps: (A) LED 1 only, (B) 
LED 2 only, and (C) LEDs 1 and 2 at the same time. The 
driver circuit consisted of a 4-bit binary counter driven 
by an Agilent 33120a pulser The Qo, Qi and Q2 out- 
puts from the counter were fed to two AND gates as 
shown in fig. |5] The logic outputs operated PNP tran- 
sistor switches that were capacitively coupled to LEDs 
1 and 2. The resulting pulses were 20 ns wide at full- 
width-half-maximum, with a variable amplitude of 0^ 
V. The pulse amplitudes, and thus the number of pho- 
tons incident on the PMT, were controlled by an Instek 
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Figure 6: Measured non-linearity of tlie PMT and electronics chain. 274 
The PMT that underwent the test was the one used in the LAB-based 275 
scintillator measurements. 

276 



GPS-3303 DC power supply, which biased the emitters 
of the two transistor switches. To reduce reflections, a 
43 Q resistor was placed in series with the LEDs, which 
had a forward bias resistance of ~6 f2. 

The voltage pulses during steps A and C, as well as 
B and C were observed to be nearly identical. The fol- 
lowing ratio R quantifies non-linearity: 



R = 100 



(P3 -fl)-f2+C 

Pi-c 



(4) 



where f 1, P2 and P3 are the ADC peak positions in 
steps A, B and C respectively, and c is the pedestal in 
number of ADC channels. The power supply voltages 
to the two transistors were varied such that P\<P2 and 
P1/P3 ~ 1/3 and P2/P3 ~ 2/3. Fig. |6] shows as a 
function of f 2 - c. The diff'erence between P3 and Pi 
is consistent with P2 at the 0.2% level, except at high 
ADC channels, when space charge effects in the PMT 
become important. The resulting gain loss makes P3 
lower than Pi + P2. We will discuss our results in terms 
of the non-linearity parameter R in the next section. 

c was measured by injecting an exponentially- 
decaying pulse of controllable amplitude into the Ortec 
571 ampHfier (fig.|2]l. It was found to be -10.72 + 0.73. 
In the process of measuring c, we confirmed the linear- 
ities of the Ortec 571 amplifier, 427 A delay and 927 
MCA units that were used to record pulse heights from 
the LS PMT. 
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Our measurements of -C{Ee) in ADC counts, for 
LAB-PPO (26 data points), Nd-doped LAB-PPO (29 
data points) and EJ-301 (18 data points), are shown 
in cross markers in figs. |7ja), (b) and (c), respectively. 
For the LAB-based scintillators, d£,/dEe, shown in the 
lower plots, increases by nearly 50 % in the energy 
range from 0.2 to 1 MeV. The change in dJ^/dE,, is less 
pronounced for EJ-301. In figs. |7ja) and (b), the circle 
markers are the 0.48, 1.04 and 4.2 MeV Compton edge 
positions (evaluated at half-height) for ys from '^^Cs, 
^"Co, and '^C(p,p')'^C, respectively. These three data 
points serve as a further verification of our alignment 
and energy reconstruction procedure. 

In ^3.21 we discussed a measurement of the non- 
linearity R (Eq. m due to the PMT and electronics at 
15 P2 values between 200 and 3500 ADC bins. This 
involved using two light sources (LEDs), one of which 
was nearly twice as strong as the other In the LS, the 
electron is also a light source. For 200<i](£'(.)<3500, we 
calculated the analogous quantity R' by interpolation: 



R' = 100 • 



[x(£.)-i:(f)]-i:(f^) + c 



(5) 



If the LS is linear with respect to E^, R' should be of 
the same order as R. The insets in figs. |7la) and (b) 
show the values of R' as a function of E^. For both LAB 
LS samples, these clearly exceed 0.2%, reaching ~20 % 
in the low-energy region. Therefore, the non-linearities 
seen in figs. |7|a) and (b) were not caused by the PMT 
and electronics. We assume this to be also the case for 
fig- Etc), although R has not been measured for the PMT 
that was used in the EJ-301 measurements. 

4.1. Interpretation and discussion 

We now provide a simple model for £.{Ee). At any 
value of Eg, the total light yield of the LS, in a wave- 
length region detectable by the PMT, is given by: 



Nt - Ns + Nc 



(6) 



where Ns is the number of photons originating from the 
scintillation process and Nc is the Cherenkov contribu- 
tion. Nc includes primary Cherenkov photons, as well 
as those emitted as a result of LS excitation by Ultra- 
Violet (UV) Cherenkov photons. JL{Ee) is proportional 
to A^T". We fitted our data to the following: 



UEe) =A-N' (£,) + B ■ Nc2iEe) + C 



(7) 



where A and B are scaling constants, and C is an offset. 
These take into account detector-related effects such as 
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Figure 7: (a) The LAB-PPO, (b) Nd-doped LAB-PPO and (c) EJ-301 electron energy scales £.(Eg) in units of ADC bins. The Compton spectrometer measurements are displayed with the cross 
markers in the upper plot. The circle markers in (a) and (b) are the Compton edge positions of y-rays from "^Cs, ^''Co, and '^C(p,p')'"C. The solid-line curves are the best fits to a simple light 
yield model (Eq.|2). The scintillation and Cherenkov components are shown as the dotted and dash-dotted curves, respectively. The lower plots show the change in dX/di?c with E^, while the 
insets show R' (defined in Eq.|3) as a function of E^. 
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photon collection and PMT quantum efficiencies, gain, 
and pedestal. The scintillation component is given by: 



r» 1 



SeiE'J 



-dEi 



(8) 



Nc2 



-ri" 



1 - 



1 



d£;dw(9) 



where fi is the electron speed in units of c, and £/ is 
defined by fi{Ef) - l/n{aj). In Eq. |9] S e is in units of 
MeV^, the frequency lj and are both in MeV, and the 
frequency integral is performed over the range 250-800 
nm. The dispersion n(w) comes from our ellipsomet- 
ric measurements of the refractive indices of SNO+ and 
EJ-301 LS above 210 nm |14]. The four fit parameters 
were A, B,C and ^g JAppendix A describes the relation 
between Eqs.|6]and|7]in more detail. 

The fit results are shown in the upper plots of fig. Q 
The scintillation and Cherenkov components are dis- 
played as the dotted and dash-dotted curves, respec- 
tively. The solid-line curves show the best fits to the 
data. For the LAB-PPO sample, the fit returned A - 
879 + 12, B = 0.78 + 0.05, C = -23.3 + 2.2, and 
kg ~ 74 //m/MeV with a x'/^^F of 12.3/22. In the case 
of Nd-doped LAB-PPO, A = 797 + 9, B = 0.77 + 0.03, 
C = -19.3 + 1.7, and kg ~ 91 yum/MeV with a x^fNT)F 
of 25.7/25. Assuming a scintillation light yield Ls of 
1 X 10"* photons/Me V, it appears that the Cherenkov light 
yield in the 250-800 nm window, Nci, had to be scaled 
up by a factor of L^BjA ~ 9 to explain our data. In 
comparison, we also fitted the data with: 



Z{E,)^A-N'{Ee) + C 



(10) 
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304 where is Birks' constant, and 5 e is the electron Stop- 350 

305 ping power, calculated using the ESTAR database [13]. 

306 Since the scintillation light yield is unknown, it was fac- 351 

307 tored into A. A^C2(£^e) was calculated with: 
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which assumes a negligible Cherenkov component, and 
a scintillation light yield governed by Birks' law. The 
X^I^Xi¥ for LAB-PPO and Nd-doped LAB-PPO were 
considerably worse: 264/23 and 492/26, respectively. 

For EJ-301, the results were A = 1558 + 9,B = 0.61± 
0.04, C = -26.9 + 3.5, and kg ~ 161 jum/MeVwith a 
X^/NDF of 14.4/14. The x^/^'DF from the fit to Eq. \T0\ 
was 51.6/15. 

Therefore, as in the KamLAND case, we find that 
the non-linearity of the SNO+ LS cannot be accounted 
for solely by Birks' law for ionization quenching. A 
combination of quenching and Cherenkov light excita- 
tion seems to better describe our measurements. For 



the SNO+ LS, the dominant contribution to the non- 
linearity appears to come from the Cherenkov compo- 
nent. This causes an increase in light yield above the 
Cherenkov threshold (-0.2 MeV). 

5. Conclusions 

To summarize, we investigated the dependence of the 
light yield on electron energy for EJ-301 and two LAB- 
based scintillators with a Compton spectrometer appa- 
ratus, whose linear response was verified. The spec- 
trometer allowed us to probe a large number of elec- 
tron kinetic energy values between 0.09 and 3 MeV. We 
found that both LAB scintillator samples demonstrate 
an enhancement in the light yield in the vicinity of the 
Cherenkov threshold for electrons. This was also ob- 
served, albeit to a lesser extent, in the EJ-301 measure- 
ments. We provided a simple model to explain these 
observations, which deviate from the usual premise of 
LS electron energy scale Unearity. 
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Appendix A. Cherenkov component 

Here we show that Nc is directly proportional to 
Nc2, which is accurately calculable with Eq. |9] Con- 
sider the two wavelength ranges: (I) A < 250 nm, 
where Cherenkov yield cannot be fully calculated be- 
cause the refractive index is unknown below 210nm, 
and (2) A>25Q nm. Nc can be decomposed as: 



Nc = fc2 ■ Nc2 + fci ■ Nci 



(A.l) 



where Nc\ and Nc2 are the numbers of Cherenkov pho- 
tons in regions 1 and 2, respectively, fci and fc2 are the 
fractions that are absorbed in regions 1 and 2, and re- 
emitted as one or more photons of longer, detectable, 
wavelengths. Since fci and fc2 have not been mea- 
sured, and Nc\ is unknown, it is not possible to perform 
an absolute calculation of Nc- However, assuming that 
the shape of the Cherenkov spectrum is independent of 
electi'on energy, ( fciNci /Nc2 + fc2) is constant at all Eg. 
Thus, Nc °^ Nc2- 



1 
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